Effect of in-situ deposition of Mg adatoms on spin relaxation in graphene 
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We have systematically introduced charged impurity scatterers in the form of Mg adsorbates to 
exfoliated single layer graphene and observe little variation of the spin relaxation times despite pro- 
nounced changes in the charge transport behavior. All measurements are performed on non-local 
graphene tunneling spin valves exposed in-situ to Mg adatoms, thus systematically introducing 
atomic-scale charged impurity scattering. While charge transport properties exhibit decreased mo- 
bility and decreased momentum scattering times, the observed spin lifetimes are not significantly 
affected indicating that charged impurity scattering is inconsequential in the present regime of spin 
relaxation times (^1 ns). 



Graphene's gate tunable transport, tabletop relativis- 
tic physics, chemical attributes, and mechanical prop- 
erties have interested researchers in a wide variety of 
fields-flH^ In particular, graphene is a candidate material 
for spintronics due to its weak hyperfine coupling and 
low intrinsic spin-orbit (SO) coupling strength {/S.so)i^^ 
which should theoretically lead to long spin lifetimes. 
Beyond scientific interest, recent progress in large area 
production by chemical vapor deposition^ combined 
with significant advances in efficient spin injection by im- 
proved tunneling contacts^ ^ has greatly improved the 
potential for advanced information processing utilizing 
spin-based logic. In particular, the introduction of effi- 
cient tunneling contacts has increased the observed spin 
lifetime by an order of magnitude (to a few ns in exfoli- 
ated graphene) by lengthening the escape time due to the 
backflow of electrons into the ferromagnetic leads. ^^E^'^ 
While graphene remains a highly promising candidate for 
carbon based spintronics, the observed spin lifetimes are 
still well below the theoretical expectations and the na- 
ture of spin relaxation remains an open question. 

In graphene, two possible spin relaxation mechanisms 
are discussed in the literatureP^^ the Elliot-Yafet (EY) 
mechanism, for which the spin relaxation time (r^) is 
proportional to the momentum scattering time (r^), and 
the D'yakonov-Perel (DP) mechanism, for which oc 
1/rp. Complicating the situation are the many possi- 
ble sources of spin relaxation in experiments on Si02 
substrate including charged impurity (CI) s catterers, 
Rashba SO coupling due toadatoms,^^^!!!^ ripples, 
and edge effects. ^^^^ Early experiments on spin trans- 
port in exfoliated graphene were able to take advan- 
tage of the tunable carrier concentration (n) and ob- 
serve a lin ear re lationship between and r^, thus sug- 
gesting EY.'^^'^ However, recent theoretical studies have 
shown that DP is expected to dominate over EY^^H^ 
and that Elliot's approach applied to graphen^^ predicts 

= {^fYtpI {^soY 1 foi" which both Fermi energy ep 
and Tp depend on carrier concentration, thus highlight- 
ing the need for experiments that can tune Tp at fixed 
n. 
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In this work we systematically introduce CI scatter- 
ers on non-local single-layer graphene (SLG) spin valves 
with high quality tunneling contacts. The experiment 
takes place in an ultra-high vacuum (UHV) deposition 
chamber with in-situ measurement capability at cryo- 
genic temperatures. All measurements and doping are 
performed in the same chamber at T=12 K and the sam- 
ple is never exposed to air. We choose Mg adsorbates as 
the CI scatterer since elements with low atomic weight 
should introduce minimal SO coupling. This substan- 
tially improves on earlier doping studies that utilized 
heavy atoms (Au) and ohmic contacts for shorter spin 
lifetimes (ts ~ 100 ps)^^, which are dominated by contact 
induced spin relaxation. We find that doping with Mg 
causes large shifts in the charge neutrality point (CNP), 
indicating significant charge transfer to the graphene 
layer, accompanied by increased momentum scattering. 
Spin transport measurements, however, indicate minimal 
effect on the spin relaxation, despite pronounced changes 
in charge transport. These results indicate that CI scat- 
tering is not an important source of spin relaxation in 
SLG in the current regime of spin lifetimes of ~1 ns. 

Graphene flakes are obtained by mechanical exfoliation 
of HOPG (SPI, ZYA) onto 300 nm SiOs/Si. SLG flakes 
are identified under an optical microscope and confirmed 
by Raman spectroscopy. The graphene flakes are electri- 
cally contacted using standard bilayer (PMMA/MMA) 
e-beam lithography and lift-off procedures. First, outer 
Au/Ti electrodes (60 nm/8 nm) are defined and de- 
posited by e-beam evaporation to serve as spin insen- 
sitive reference contacts. The sample is then annealed 
for 3 hours in UHV at 150° C immediately prior to the 
second lithography step, which defines the inner ferro- 
magnetic electrodes. Angle evaporation is utilized to de- 
posit sub-monolayer Ti02, which serves as a diffusion 
barrier for the 0.9 nm MgO tunnel barrier, and 80 nm 
Co. These tunneling contacts are deposited in a molecu- 
lar beam epitaxy (MBE) chamber with base pressure of 
1 X 10~^^ torr. The electrodes are then capped with 5 
nm AI2O3. A detailed description for the fabrication of 
tunneling contacts is described elsewhere. 

Charge and spin transport measurements at T=12 K 
are performed on non-local devices as shown in Fig. [l^. 
The gate dependent resistivity of pristine SLG (sam- 
ple A) is shown in Fig. lib with maximum resistivity at 
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FIG. 1. (a), Device schematic of the non-local spin valve 
geometry with inner Co electrodes (blue) and outer Au elec- 
trodes (yellow), (b) Gate dependent resistivity for sample 
A at T = 12 K. (c) Rnl for pristine SLG at Vg = V. 
(d) Hanle spin precession data in parallel (black) and anti- 
parallel (red/grey) configuration between electrodes B and C 
for pristine SLG at Vb = V. A constant spin-independent 
background has been subtracted. 



the charge neutrality point, Vcnp = —20 V. The mo- 
bility is calculated by taking the slope of the conduc- 
tivity (/i = Aa/eAn) where the carrier concentration, 
n (positive for holes), is determined using the relation 
n = -a{VG - Vcnp) and a = 7.2 x 10^° V-^cm-^ for 
300 nm Si02 gate dielectric. The resulting electron and 
hole mobilities are /ig = 1774 cm^/Vs and jHh = 1508 
cm^/Vs, respectively. For spin transport measurements, 
an AC current, Inl = ^ A^A (11 Hz), is applied to inject 
spin-polarized carriers into SLG at electrode B. This spin 
polarization diffuses through the graphene channel along 
the X-axis to electrode C. A non-local voltage, Vnl^ is 
detected using standard lock-in techniques between elec- 
trodes C and D due to the accumulation of spins be- 
neath electrode C. The detected voltage, Vatl, is propor- 
tional to the spin-dependent chemical potential difference 
between electrodes C and D.^ The non-local resistance, 
Rnl = Vnl/Inl^ depends on the relative orientation 
of the two inner ferromagnetic electrodes and is posi- 
tive (negative) for parallel (antiparallel) alignment. An 
external magnetic field, Bapp^y, is applied along the elec- 
trode easy axis (?/-axis) and is used to control the relative 
alignment of the magnetic electrodes. A typical sweep of 
Bapp,y for sample A at Vb = V (n = -1.44 x 10^^ 
cm~^) is shown in Figure [l]^, for which the spin signal 
ARnl = Rnl ~ Rnl i^ ^^-^ ^- dimensions of the 
graphene spin channel for sample A are defined by the 
channel length L = 2.2 /am and width w = 2A /im. The 
spin lifetime can be determined from Hanle spin preces- 



sion measurements in which a magnetic field, Bapp^z-, is 
applied out of plane allowing the injected spins to precess 
around Bapp^z- At large fields, the ensemble spin popula- 
tion dephases as Bapp^z is increased due to a distribution 
of arrival times at electrode C. In the tunneling limit, 
the ensemble spin precession can be fit using the Hanle 
equation,'^^ 

Rnl oc / — cos ^ t ] e ' ""at (1) 

7o V47rDt V ^ / 

where D is the diffusion coefficient, g is the electron 
^-factor, jiB is the Bohr magneton, and h is the re- 
duced Planck's constant. Figure [T]i shows characteristic 
Hanle curves for parallel and antiparallel alignment for 
n = — 1.44 X 10^^ cm~^, where best fits to the Hanle equa- 
tion yield the diffusion coefficient D = 0.058 m^/s, spin 
lifetime = 1.10 ns, and spin diffusion length = 8.0 
/im. 

Next, Mg adsorbates are deposited in- situ in the UHV 
MBE chamber with base pressure 3 x 10~^^ torr while the 
sample is maintained at T=12 K. Elemental Mg (99.99%) 
is evaporated from an effusion cell at a rate of 0.055 
A/min calibrated by a quartz crystal monitor and corre- 
sponds to a doping rate of 0.02% of a monolayer (ML) per 
second, where 1 ML is defined as 1.908 x 10^^ atoms/cm^. 
After 1 s Mg deposition, the charge and spin transport 
are re-measured. Figure [2] summarizes the effect on the 
charge transport on sample A following Mg doping. Fig- 
ure [2^ shows conductivity a vs. Vq for Mg doping of 
sample A up to 7 s deposition time. After 7 s of Mg dop- 
ing, Vcnp has shifted to Vg = — 70 V. This indicates that 
Mg donates electrons to the graphene, consistent with re- 
ports on transition metals and potassium. ^^ ^^ Figure [2]3 
displays Vcnp for each doping time and demonstrates a 
linear relation between charge transfer and Mg coverage 
at a rate of —1438 V/ML. Also, Mg doping introduces 
CI scattering which decreases the conductivity and the 
mobility. Figure [2]3 displays the effect of systematic Mg 
doping on the electron and hole mobilities. For undoped 
graphene, the mobility is ji^ = 1774 cm^/Vs and and 
jih = 1508 cm^/Vs, and decreases to /ig = 599 cm^/Vs 
and fih = 453cm^/Vs after 7 s deposition time. The mo- 
mentum scattering time can be determined using Boltz- 
mann transport theory,'^ 
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where h is Planck's constant, e is the electron charge, 
'^i? ~ 1 X 10^ m/s is the Fermi velocity, and ge = 2 and 
gy = 2 are the spin and valley degeneracies. Fig. [2]i 
shows Tp vs. Mg doping for electrons and holes at n = 
±1.44 X 10^^ cm~^. With increasing Mg coverage, the 
momentum scattering time decreases due to increased 
CI scattering. Lastly, we investigate the nature of Mg 
morphology on the graphene surface. Figure [2^ shows 
the shift in Dirac point plotted against l/fiavg — I/mo, 
where jiavg is the average of the electron and hole mo- 
bilities and /io is the average electron and hole mobility 
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FIG. 3. (a) Hanle spin lifetime for sample A plotted against 
Mg coverage for electrons (red/grey) and holes (black). Inset 
is the diffusion coefficient obtained from the Hanle fit as a 
function of Mg coverage at T=12 K. (b) Hanle spin lifetime for 
sample A plotted against the calculated momentum scattering 
times for each doping for both electrons (red/grey) and holes 
(black), (c) Electron (blue/grey) and hole (black) mobility 
as a function of Mg coverage for sample B at T=12 K. (d) 
Hanle spin lifetime for sample B plotted against the calculated 
momentum scattering times. 



FIG. 2. (a) Gate dependent conductivity as the SLG (sample 
A) is systematically doped up to 7 s with Mg adsorbates at 
T=12 K. (b) Charge neutrality point (CNP) plotted against 
the Mg doping time, (c) Electron (red/grey) and hole (black) 
mobility as a function of Mg coverage, (d) Calculated mo- 
mentum scattering time for electrons (red/grey) and holes 
(black) as a function of Mg coverage, (e) Shift in the charge 
neutrality point, — AVcatp, plotted against the change in in- 
verse mobility. The dashed line is a power law fit (best fit 
exponent b = 0.72). 



for pristine graphene. The dashed line is a power law 
fit of —AVcNP {l/fiavg — 1/Mo)^7 ^^t which values 
of 1.2 < b < 1.3 indicates a 1/r scattering potential for 
point-like scatterers.^^ The best fit value of 6 = 0.72 
suggests the possibility of clustering even at cryogenic 
temperatures. ^^-^^ This does not introduce a theoretical 
difficulty because the relationship = {epYTp/ {/S^soY 
for EY scattering in SLG has been shown to hold for a 
wide variety of scattering sources including weak scat- 
terers, strong scatterers (i.e. vacancies), CI scatterers, 
and clusters.'^ Lastly, we note that the gate dependent 
resistance curves exhibited no measurable change as a 
function of time in between Mg depositions. 

We now turn to the effect on spin relaxation in SLG 
by Mg doping. After each Mg deposition at 1 s intervals, 
Hanle spin precession measurements were performed for 
n = ±1.44 X 10^^ cm~^. The resulting fits to the Hanle 
curves yield values for and D which are plotted against 



Mg doping time in Figure [3^ and|3^ inset, respectively. 
The diffusion coefficient decreases with increasing Mg 
coverage starting at 0.058 m^/s (0.038 m^/s) for pristine 
graphene and decreases to 0.021 m^/s (0.033 m^/s) for 7 
s doping time for electrons (holes). This is in agreement 
with the observed charge transport behavior for which 
momentum scattering increases with Mg doping. Inter- 
estingly, the spin lifetimes (Fig. [3|i) show minimal varia- 
tion, without a significant trend for electrons and holes. 
In Figure we plot the Hanle spin lifetime for sam- 
ple A against the momentum scattering time calculated 
from the conductivity using equation Q from Boltzmann 
transport theory for sample A. With increasing Mg dop- 
ing, Tp decreases from ~ 35 fs to ~ 20 fs, but the spin 
relaxation time is constant for holes (black squares) while 
decreasing only slightly for electrons (red/grey squares). 
This experiment was repeated on several samples and 
in general does not display any substantial variations 
as a function of r^. For instance, results for a sample 
with higher initial mobility (sample B) are summarized in 
Fig.jSj) and 3d. Figure [Sj: displays the change in mobility 
for electrons and holes under Mg doping. For sample B, 
the mobility decreases from 4415 cm^/Vs (3200 cm^/Vs) 
for the pristine spin valve to 598 cm^/Vs (1290 cm^/Vs) 
after 7 s Mg doping for electrons (holes). In Fig. [sji, 
we show Ts displayed against the momentum scattering 
times for sample B at n = ±1.44 x 10^^ cm~^. Here, 
Ts is near 800 ps and stays relatively unchanged as Tp 
decreases from ~ 60 fs to ~ 24 fs. 
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As Fig. [3]3 and 3d show, does not vary substantially 
as Tp is varied by CI scattering. This is in agreement 
with and goes beyond recent reports on CI scattering by 
organic-ligand bound nanoparticles, which are able to re- 
versibly tune the mobility and momentum scattering.^ 
Due to the relatively large size (~ 13 nm, which is over 50 
lattice constants) of the nanoparticles used in that study, 
it is not possible to draw conclusions for atomic-scale CI 
scatterers such as surface adatoms and impurities in the 
Si02 substrate. In contrast, Mg adsorbates are able to 
probe the atomic-scale regime. With the agreement be- 
tween two quite different experiments (Mg adsorbates de- 
posited in UHV and organic-ligand bound nanoparticles 
deposited by drop casting) probing different length scales 
of the scattering potential, it is clear that spin relaxation 
in graphene is not determined by CI scattering despite 
its importance for momentum scattering. 

It is also worth mentioning that this result is not in- 
compatible with the early experiments showing a lin- 
ear relation between and D by tuning the carrier 
concentration.^^ While Mg adsorbates modify Tp by 
introducing CI scattering and possibly local Rashba SO 
coupling, there are many alternative sources which might 
contribute to EY (i.e. weak scatterers, resonant scatter- 
ing, phonon scattering) which could still present them- 
selves as the carrier concentration is tuned leading to 



Tg ^ D. Thus, EY spin relaxation originating from 
sources other than CI scattering is still viable. 

Some other possibilities to consider are DP spin relax- 
ation in spatially inhomogeneous Rashba SO fields. It 
has recently been proposed that this type of SO coupling 
can result in a competition between EY-like and DP-like 
scaling behavior to yield unconventional scaling between 
Tg and Tp}^ Another possibility is that the spin lifetime 
is limited by contact effects such as inhomogeneous stray 
fields.!^ Due to its atomically thin nature, this could have 
a larger effect for graphene compared to semiconductor or 
metallic spin transport systems that are typically much 
thicker. 

In conclusion, we have investigated charge and spin 
transport in SLG by systematically introducing CI scat- 
terers in the form of Mg adsorbates. The introduction 
of Mg was shown to transfer electrons to the SLG and 
decrease the momentum scattering time. Despite pro- 
nounced changes in momentum scattering, no significant 
variation was seen in spin relaxation. This indicates that 
CI scattering is not an important source of spin relax- 
ation in SLG in the current regime of spin lifetimes (~1 
ns). 
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NRI-NSF (NEB-1124601). 



^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. L Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 

A. A. Firsov, Nature 438, 197 (2005). 

^ Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 

438, 201 (2005). 
^ A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009). 

^ C. Lee, X. Wei, J. W. Kysar, and J. Hone, Science 321, 
385 (2008). 

^ D. Huertas-Hernando, F. Guinea, and A. Brataas, Phys. 

Rev. B 74, 155426 (2006). 
^ H. Min, J. E. Hill, N. A. Sinitsyn, B. R. Sahu, L. Kleinman, 

and A. H. MacDonald, Phys. Rev. B 74, 165310 (2006). 
^ Y. Yao, F. Ye, X. L. Qi, S. C. Zhang, and Z. Fang, Phys. 

Rev. B 75, 041401 (R) (2007). 
^ N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and 

B. J. van Wees, Nature 448, 571 (2007). 

^ X. Li, W. Cai, J. An, S. Kim, J. Nah, D. Yang, R. Finer, 
A. Velamakanni, I. Jung, E. Tutuc, S. K. Banerjee, 
L. Colombo, and R. S. Ruoff, Science 324, 1312 (2009). 
^° S. Bae, H. Kim, Y. Lee, X. Xu, J.-S. Park, Y. Zheng, 
J. Balakrishnan, T. Lei, H. R. Kim, Y. I. Song, Y.-J. Kim, 
K. S. Kim, B. Ozyilmaz, J.-H. Ahn, B. H. Hong, and 
S. Ijima, Nature Nano. 5, 574 (2010). 

W. Han, K. Pi, K. M. McCreary, Y. Li, J. J. I. Wong, 

A. G. Swartz, and R. K. Kawakami, Phys. Rev. Lett. 
105, 167202 (2010). 

B. Dlubak, M. B. Martin, C. Deranlot, B. Servet, S. Xavier, 
R. Mattana, M. Sprinkle, C. Berger, W. A. D. Heer, 
D. Petroff, A. Anane, P. Seneor, and A. Fert, Nature Phys. 
8, 557 (2012). 



^2 H. Dery, H. Wu, B. Ciftcioglu, M. Huang, Y. Song, R. K. 
Kawakami, J. Shi, I. Krivorotov, I. Zutic, and L. J. Sham, 
IEEE Trans. Electron. Dev. 59, 259 (2012). 

14 ^ _Y. Yang, J. Balakrishnan, F. Volmer, A. Avsar, 
M. Jaiswal, J. Samm, S. R. Ah, A. Pachoud, M. Zeng, 
M. Popinciuc, G. Giintherodt, B. Beschoten, and 

B. Ozyilmaz, Phys. Rev. Lett. 107, 047206 (2011). 

W. Han and R. K. Kawakami, Phys. Rev. Lett. 107, 
047207 (2011). 

C. Ertler, S. Konschuh, M. Gmitra, and J. Fabian, Phys. 
Rev. B 80, 041405(R) (2009). 

F. Ochoa, A. H. Castro Neto, and F. Guinea, Phys. Rev. 
Lett. 108, 206808 (2012). 

A. H. Castro Neto and F. Guinea, Phys. Rev. Lett. 103, 
026804 (2009). 

P. Zhang and M. W. Wu, Phys. Rev. B 84, 045304 (2011). 

D. Huertas-Hernando, F. Guinea, and A. Brataas, Eur. 
Phys. J. Spec. Top. 148, 177 (2007). 

D. Huertas-Hernando, F. Guinea, and A. Brataas, Phys. 
Rev. Lett. 103, 146801 (2009). 

C. Jozsa, T. Maassen, M. Popinciuc, P. J. Zomer, 
A. Veligura, H. T. Jonkman, and B. J. van Wees, Phys. 
Rev. B 80, 241403(R) (2009). 

M. Popinciuc, C. Jozsa, P. J. Zomer, A. Veligura, H. T. 
Jonkman, and B. J. van Wees, Phys. Rev. B 80, 214427 
(2009). 

2^ P. Zhang and M. W. Wu, New J. Phys. 14, 033015 (2012). 

C. Weeks, J. Hu, J. Ahcea, M. Franz, and R. Wu, Phys. 

Rev. X 1, 021001 (2011). 
2^ K. Pi, W. Han, K. M. McCreary, A. G. Swartz, Y. Li, and 

R. K. Kawakami, Phys. Rev. Lett. 104, 187201 (2010). 

G. Giovannetti, P. A. Khomyakov, G. Brocks, V. M. 



5 



Karpan, J. van den Brink, and P. J. Kelly, Phys. Rev. 
Lett. 101, 026803 (2008). 

K. Pi, K. M. McCreary, W. Bao, W. Han, Y. F. Chiang, 
Y. Li, S. W. Tsai, C. N. Lau, and R. K. Kawakami, Phys. 
Rev. B 80, 075406 (2009). 

K. M. McCreary, K. Pi, A. G. Swartz, W. Han, W. Bao, 
C. N. Lau, F. Guinea, M. L Katsnelson, and R. K. 
Kawakami, Phys. Rev. B 81, 115453 (2010). 
J. H. Chen, C. Jang, S. Adam, M. S. Fuhrer, E. D. 
Williams, and M. Ishigami, Nature Phys. 4, 377 (2008). 
Y. W. Tan, Y. Zhang, K. Bolotin, Y. Zhao, S. Adam, E. H. 



Hwang, S. D. Sarma, H. L. Stormer, and P. Kim, Phys. 
Rev. Lett. 99, 246803 (2007). 

S. Adam, E. H. Hwang, V. M. Galitski, and S. D. Sarma, 
PNAS 104, 18392 (2007). 

W. Han, J.-R. Chen, D. Wang, K. M. McCreary, H. Wen, 
A. G. Swartz, J. Shi, and R. K. Kawakami, Nano Lett. 
12, 3443 (2012). 

S. P. Dash, S. Sharma, J. C. L. Breton, J. Peiro, H. Jaffres, 
J.-M. George, A. Lemaitre, and R. Jansen, Phys. Rev. B 
84, 054410 (2011). 



